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Abstract: A rich variety of brane cosmologies is obtained once one allows for energy 
exchange between the brane and the bulk, depending on the precise form of energy trans- 
fer, on the equation of state of matter on the brane and on the spatial topology. This 
is demonstrated in the context of a non-factorizable background geometry with zero ef- 
fective cosmological constant on the brane. An accelerating era is generically a feature 
of these solutions. In the case of low-density flat universe more dark matter than in the 
conventional FRW picture is predicted, while spatially compact solutions are found to 
delay their re-collapse. In addition to the above, which the interested reader will find in 
greater detail in (!]], a first attempt towards a complete description of the full dynamics 
of both the bulk and the brane is reported. 



1. Introduction 

The idea that we might be living inside a defect, embedded in a higher dimensional space 
has already a long history. Concerning the nature of the defect, a solitonic codimension 
one or higher topological object was proposed [§] in the context of an ordinary higher 
dimensional gauge field theory, coupled ||] or not to gravity. It was soon realized, however, 
that, in contrast to scalar and spin-1/2 fields, it would be difficult to confine gauge fields on 
such an object. Various interesting ideas and scenaria were studied [||, which have not yet 
given a fully satisfactory picture. In connection with the topology and the size of the bulk 
space on the other hand, the popular choice was that the extra dimensions are compact, 
with size of order O(Mp^). It was argued in ||, that in the context of the heterotic string 
with supersymmetry broken a la Scherk-Schwarz this should not be true anymore. The 
scale of supersymmetry breaking is tied to the size of internal dimensions, and a desirable 
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supersymmetry breaking scale of a few TeV implies an extra dimension of about 10 -16 cm. 
Despite difficulties to build a realistic model based on these ideas, the scenario was taken 
seriously and analyzed further for its phenomenological consequences ||. 

The situation is drastically different in the context of type-I string theory |7j. A few 
developments led to an exciting possibility and renewed interest in the whole idea. First, 
with the discovery of the D-branes as an essential part of the "spectrum" in type-I string 
theory, one could conjecture that we inhabit such a D-brane embedded in a ten-dimensional 
bulk. The usual solitonic defect of field theory was thus replaced by an appropriate col- 
lection of D-branes, which by construction confine the gauge fields ||, together with all 
the ingredients of the standard model. All known matter and forces lie on our brane 
world [k|, with the exception of gravity, which acts in the bulk space as well. It was, 
however, pointed out [11], that for Kaluza-Klein extra dimensions the gravitational force 
on the brane was consistent with all laboratory and astrophysical experimental data, as 
long as the extra dimensions were smaller than a characteristic scale. This led to the ex- 
citing possibility of two extra dimensions in the sub-millimeter range. Furthermore, it was 
demonstrated in the context of an appropriate effective five-dimensional theory of gravity, 
that once we take into account the back reaction of the brane energy-momentum onto the 
geometry of space-time, the graviton is effectively confined on the brane and Newton's law 
is reproduced to an excellent accuracy at large distances, even with a non-compact extra 
dimension [12|. 

At the same time, the analysis of the cosmological consequences of the above hypothe- 
ses attracted considerable interest. The first step was taken in jDJ, where the evolution of 
perfect fluid matter on the brane was studied, with no reference to the bulk dynamics and, 
therefore, no energy transfer between the brane and the bulk. Alternatively, a bulk-based 
point of view was adopted in 14], where the cosmology induced on a moving 3-brane in 
a static Schwarzschild-AdSs background was studied, and a general interpretation of cos- 
mology on moving branes, together with the idea of "mirage" cosmology, were presented 
in |ll|. The equivalence of the two approaches was demonstrated in ||16|| . It was also 
pointed out that branes provide natural mechanisms for a varying speed of light [17]. 

Energy-exchange between the brane and the bulk should in principle be included in any 
realistic cosmological scenario, and its effects have been studied in detail in the context of 
flat compact extra dimensions [18]. The role of energy-exchange on brane cosmology in the 
case of non-factorizable extra dimensions has not yet been investigated extensively [19|, 20], 
even though the importance of energy outflow from the brane has already been demon- 
strated in the context of a Randall- Sundrum configuration with additional gravity induced 
on the brane pC[[. 

Here an attempt towards a more complete analysis of the cosmic evolution of the brane 
in the presence of energy flow into or from the bulk is presented. The aim is to generalize 
the picture in bulk AdS (bulk gravity plus cosmological constant) by considering a general 
bulk theory (that includes gravity). There may be more bulk fields and more general bulk- 
brane couplings. As argued below, one may consider a regime where the bulk energy can 
be consistently neglected from the equations. Moreover, the energy exchange component 
of the energy-momentum tensor, will be assumed to be just a power of the matter density 
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of the brane. Although there are other valid parameterizations, the one chosen here is 



motivated by our previous work in [ 20 1 , where we analyzed the out-flow of energy from the 
brane due to graviton radiation in the presence of induced gravity. Standard cross section 
calculations give an outflow that is a function of the temperature and other fundamental 
constants of the theory. If one re-expresses the temperature in terms of the running density 
using the cosmological equations we obtain a rate of flow that is a power of the density, 
with a dimensionful coefficient that depends on fundamental constants and initial energy 
densities. This is valid typically for a whole era, that is, piecewise in the cosmological 
evolution. 

In a given theory, with a specific bulk content and brane-bulk couplings, the exponent 
of the density as well as the coefficient are calculable functions of the coupling constants 
of the Lagrangian as well as initial densities. It is also obvious that if the bulk theory 
is approximately conformal the most general form of rate of outflow will be polynomial 
in the density. In any case, the analysis carried out here is " phenomenological" , and the 
parameters used in the energy exchange function have not been derived from fundamental 
dynamics. 

The presentation is organized in five sections of which this introduction is the first. 
In section [2] the framework of this talk is described and the approximations on the brane- 
bulk exchange are presented. The effective set of equations for the analysis of the brane 
cosmology are derived. Section ^ contains several interesting characteristic solutions of 
the brane cosmology. In Section |3] a few brief comments will be offered concerning the 
exact treatment of the influx/outflow equations, together with a couple of representative 
cosmological histories, plotted to show the rich variety one may expect in the context of 
this scenario. Section ^ contains a first attempt to describe the full dynamics of both 
the bulk and the brane. With appropriate perfect fluids in the bulk and on the brane 
one may describe the late time cosmological evolution of the brane-world. Our results are 
summarized and the prospects for further research along these lines are discussed in the 
final section. 

2. The model 

The basic assumption is that we live on a 3-brane embedded in a 4+1-dimensional bulk 
with the transverse dimension assumed non-compact. In particular, we shall be interested 
in the generic model described by the action 

S = Jd 5 x {M 3 R - A + £g at ) + J rf 4 x^ [-V + ££ iat ) , (2-1) 

where R is the curvature scalar of the five-dimensional metric qab, A,B = 0, 1, 2, 3, 5, A is 
the bulk cosmological constant, and g a j3i with a, /3 = 0, 1, 2, 3, is the induced metric on the 
3-brane. We identify (x, z) with (x, —z), where z = x§. However, following the conventions 



of [12] we extend the bulk integration over the entire interval (— oo,oo). The quantity 
V includes the brane tension as well as quantum contributions to the four-dimensional 
cosmological constant. £™ at represents the matter on the brane. It includes massless 
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excitations (branons) related to fluctuations of the brane, which is assumed frozen at z = 0. 
Finally, £g at stands for the bulk matter action. 

As explained in the introduction, various instances of this model have been extensively 
discussed in the literature. Here we shall present some aspects of the influence of energy 
exchange between the brane and the bulk upon the cosmological evolution of our world on 
the brane. 

We consider an ansatz for the metric of the form 

ds 2 = -n 2 (t, z)dt 2 + a 2 {t, z)- fij dx i dx j + b 2 (t, z)dz 2 , (2.2) 

where jij is a maximally symmetric 3-dimensional metric. We use k = —1,0, 1 to parame- 
terize the spatial curvature. 

The non-zero components of the five-dimensional Einstein tensor are 

fd/d b\ n 2 fa" a' fa' b'\\ 7 n 2 \ 
'' '-' M a I a n j b In a j a n J 



b 2 

a 2 \ a ( d „n\ 'a b ( 'a h\ b \ , , n A . 

+ ^{- a {-a +2 n)- 2 a + l K + n) ~ l) ~ ^ 

a.^Jti + tk fL), (2 .5) 

\ n a a b a J 

^_ 3 {*(s^)_4(*(*_*) + *)_ t 4}. (2 . 6) 

l a \ a n J n z \a \a n J a J or ) 

Primes indicate derivatives with respect to z, while dots derivatives with respect to t. 
The five-dimensional Einstein equations take the usual form 

GAC= 2^3 TAC ' (2 - 7) 

where Tac denotes the total energy- momentum tensor. 

Assuming a perfect fluid on the brane and, possibly an additional energy-momentum 
Tq \m,B in the bulk, we write 

T C = Tcl, b + T A c \ m>b + T A c \ vB + T A c \ mB (2.8) 

T cl, b = ^diag(-v,-v;-v;-v;o) (2.9) 



T A d tB = diag(-A, -A, -A, -A, -A) (2.10) 
T cLb = ^diag(-p,p,p,p,0), (2.11) 



where p and p are the energy density and pressure on the brane, respectively. The behavior 
of T^\ m ^B is in general complicated in the presence of flows, but we do not have to specify 
it further at this point. 
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We wish to solve the Einstein equations at the location of the brane. We indicate by 
the subscript o the value of various quantities on the brane. Integrating equations (|2.3|), 



fl2.4j ) with respect to z around z = gives the known jump conditions 



i + = = -J^p b o a o (V + p) (2.12) 

i+ = -n' - = -^m h ° n ° { ~ v +2p+ 3p) ■ (2 - 13) 



The other two Einstein equations (|2.5[) and (|2.6|) give 



n' a a' b a' 1 
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A ^o + 7777^55 , (2-15) 



N3 

where ^5^55 are the 05 and 55 components of TAc\m,B evaluated on the brane. Substi- o 
tuting ( |2~lp , ( gig ) in equations ( g3g ), d2~T5|) one obtains 

■ 

p + 3^(p + p) = -fan (2.16) 

2. (ao foo] 2 _ aonA k_ J_ / ^ y2 \ _ (2 u) 

n 2 \a \a J a n )^a? QM*\ ) { ' 



We are interested in a model that reduces to the Randall-Sundrum vacuum [12] in the 
absence of matter. In this case, the first term on the right hand side of equation ( 2.17| ) 
vanishes. A new scale kns is defined through the relations V = —A/kjis = 12M 3 £;rs. 

We shall further assume that the energy exchange between the bulk and the brane has 
negligible effect on the bulk. With this assumption, which is equivalent to assuming that 
the brane is embedded in a "heat reservoir", we shall be able to derive a solution of the 
brane cosmology that is independent fo the bulk dynamics. 

At this point we find it convenient to employ a coordinate frame in which b Q = n = 1 
in the above equations. This can be achieved by using Gauss normal coordinates with 
b(t, z) = 1, and by going to the temporal gauge on the brane with n = 1. The assumptions 
for the form of the energy-momentum tensor are then specific to this frame. 1 



1 If the vacuum energy dominates over the matter content of the bulk, we expect that the form of the 
metric will be close to the Randall-Sundrum solution with a static bulk. Thus, we expect (even though we 
cannot demonstrate explicitly without a full solution in the bulk) that in a generic frame, in which 



rrnldiag 

1 lm,B 



rrildiag 

1 \v,B 



< 1 (2.18) 



we shall have b ~ 0. Then the transformation that sets fo = 1 is not expected to modify significantly the 
energy-momentum tensor. 



- 5 - 
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Using f3 = M 6 /144 and 7 = VM 6 /144, and omitting the subscript o for convenience 
in the following, we rewrite equations ( 2.16| ) and ( |2.17| ) in the equivalent form 



p + 3(1 + w)-p = —T (2.19) 
a 

^ = PP 2 + 27P - 4 + X ( 2 - 2 0) 



X + ^ x = 2(i(p+l)T, (2.21) 

where p = wp, T = 2T° 5 is the discontinuity of the zero-five component of the bulk energy- 
momentum tensor. As mentioned above, with an appropriate choice of parameters we have 
set to zero the effective cosmological constant A = (A + U 2 /12M 3 )/12M 3 on the brane, 
which should otherwise be added in th eright-hand side of ( |2.20| ). 



Equations Q2.1S| ), ( |2.20D and ( 2.21 ), in place of the usual FRW ones, supplemented by 
a set of initial conditions (pi,a>i,Xi) or, equivalently, (pi,a,i,ai) 2 , describe the dynamical 
evolution of matter on the brane, within the assumptions stated above. They have a 
straightforward interpretation. The energy density evolution equation ( 2.1S| ) is modified 



by the presence of the energy exchange term in the right-hand side. Eq. ( |2.20| ) is the 
modification of Friedman equation. It has the extra term quadratic in the energy density, 
and may be thought of as the definition of the auxiliary density x- With this definition 
the other two equations are equivalent to the original system ( 2.16| ), ( [2.17 ). As we will 



see later on, in the special case of no-exchange (T = 0) % represents the mirage radiation 
reflecting the non-zero Weyl tensor of the bulk. 

The second order equation ( p. 17 ) for the scale factor becomes 



- = -(2 + 3w)f3p 2 - {l + 3w)jp- X - (2.22) 
a 

Notice that in the special case of iy = 1/3 one may define a new function x = X + 27/3. 
The functions x, p and a satisfy equations ( |2.19| ) to ( |2.22| ) with x 111 place of x an d 7 = 0. 
This should be expected, since for w = 1/3 there is no 7 left in equation ([2.170 . 

It is not surprising that the introduction of the bulk and its interaction with the brane 
enriches the possible cosmologies on the latter. The interpretation of observations on the 
brane becomes unavoidably more complicated, as a consequence of the at best indirect 
knowledge of the evolution of the bulk and of its interaction with brane matter. This will 
become clear in the discussion of possible solutions to which we turn next. 

3. Special solutions 



Before we embark on the discussion of the general solutions of equations fl2.19| )-( ^.21 ), it 



is instructive to start with a few special cases easy to treat analytically and whose physical 



2 In the present case {pi, at) of conventional FRW cosmology are not enough. The reason is that the 



generalized Friedmann equation (2.2C) (or (3.5)) is not a first integral of the Einstein equations because of 



the possible energy exchange between the brane and the bulk. 
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content is more transparent. The analysis simplifies considerably in the low-density region, 
defined by p <C 7/ (3 = V , the case in which the matter energy density on the brane is much 
smaller than the brane cosmological term. In this case one may ignore the term (3p 2 in the 
above equations compared to 7/3. As a result, equations ( |2.19 )-( 2.21 ) reduce to 



p + 3(1 + w)Hp = -T (3.1) 

H 2 ={-) 2 = 2 7 p + x-4 M 
\a J gt 

X + 4Hx = 2 7 T. (3.3) 



3.1 Constant acceleration with w 7^ — 1 

An interesting feature of this framework is the possible presence of accelerating cosmological 
solutions with constant acceleration, even for w/-l, that would be required in the context 
of the standard FRW cosmology. We can have exponential expansion with constant Hubble 
parameter H = H* , even if the brane content is not pure vacuum energy. 



Indeed, let us choose for simplicity k = and search for a solution of equations (3.1)- 
Q3.3p with constant H = H*. It is straightforward to show that in this case the density p 
and the auxiliary field % are a ^ so constant 3 . 

Denote with a * subscript the constant values of these quantities. They satisfy 

3H(l + w)p m = -T(p m ) (3.4) 
Hi = 2 1P * + X * (3.5) 
2H* X * = lT(p*). (3.6) 

It is clear from equation (|3.4D that, for positive matter density on the brane (p > 0), 
flow of energy into the brane (T(p) < 0) is necessary. The accretion of energy from the 
bulk depends on the dynamical mechanism that localizes particles on the brane. It is not 
difficult to imagine scenaria that would lead to accretion. If the brane initially has very low 
energy density, energy can by transferred onto it by bulk particles such as gravitons. An 
equilibrium is expected to set in if the brane energy density reaches a limiting value. As a 
result, a physically motivated behavior for the function T(p) is to be negative for small p 
and cross zero towards positive values for larger densities. In the case of accretion it is also 
natural to expect that the energy transfer approaches a negative constant value for p — > 0. 
The solution of equations ( |3.4[) -( |3^) satisfies 

Tfo,) = -^(1 + w)(l- 3wy/ 2 p 3 J 2 (3.7) 

Hi = (3.8) 
3(1 + w) 

X* = o 7P* • (3-9) 



3 Combining (3.1) and (3.3) one obtains \ + + 67(1 + w)H*p + AH,\ 
(3.2) gives x = — 2p/(37(l + w)). This, upon substitution into (3.2) leads to p 



= 0, which, combined with 
^constant. 
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What is the interpretation of this result? For a general form of T(p), determined by the 
physics of fundamental interactions, equation is an algebraic equation with a discrete 
number of roots. For any value of uu in the region —1 < w < 1/3 a solution is possible. The 
resulting cosmological model is characterized by a scale factor that grows exponentially 
with time. The energy density on the brane remains constant due to the energy influx 
from the bulk. The model is very similar to the steady state model of cosmology [?]. 
The main differences are that the energy density is not spontaneously generated, and 
the Hubble parameter receives an additional contribution from the "mirage" field x ( see 
equation (|3.5|)). 

By linearizing equations fl3.1| )-([0|) around the above solution, one may study its stabil- 
ity. The conclusion is that for — 1 < w < 1/3 and < v < 3/2, where v = din \T\(p*)/dln p, 
the fixed point solution described here is stable 

For w = — 1 we get the standard inflation only for a value that is a zero of T{p). In 
this case there is no flow along the fifth dimension and also x* = 0- 



3.2 "Mirage" radiation for energy outflow 

Another interesting situation arises in the case of outflow T > 0, of energy from the 
brane into the bulk and with radiation domination on the brane, i.e. with p = p/3 4 . 
Equations ( |3.1[ ) and ( |3.3| ) have an exact solution independent of the explicit form of T: 

' + £ = («+!) i < 3 - 10 > 

and 

H 2 = {2 1Pi + x.) °\ - 4 • (3-11) 

Assume that initially xi = 0- It is clear that the effect of the radiation on the expansion does 
not disappear even if it decays during the cosmological evolution: the Hubble parameter 
of equation ( 3.1 1| ) is determined by the initial value of the energy density, diluted by the 



expansion in a radiation dominated universe. The real radiation energy density p, however, 
falls with time faster than a~ 4 . 

To appreciate the unusual features of the resulting cosmology, consider the simple case 
of what we shall call "radioactive" matter on the brane, T = Ap with A > 0. Then 



equation (3.1) can be integrated for arbitrary w to give 



e~ M , (3.12) 

where we have considered the general case p = wp. This equation has an obvious interpre- 
tation. The energy density on the brane dilutes both as a result of the expansion of the 
universe (the a~ 3 ( 1+w ^ factor) as well as a consequence of its "radioactive decay" (the e~ At 



4 This is the simplest case to treat analytically. The conclusions of this subsection generalize Q| to other 
interesting cases, such as dust on the brane. 
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factor). In the special case ro = 1/3 we are considering here, one obtains 



a i -At 

ar 



(3.13) 



X = 2 1Pi ^ (1 



-At 



(3.14) 



Imagine an initial energy density pi ^ 0, while x-i = 0- The Hubble parameter, given by 
equation ( 3.1l| ) with Xi = 0j corresponds to an initial radiation density pi, further diluted 
only by the expansion. At late times the energy density disappears exponentially fast and 
the expansion is the consequence of a "mirage" effect, sustained, through the auxiliary \i 
by the expansion itself. Of course, the presence of a "mirage" term is possible even without 
energy flow, since equation ( |3.3| ) has a solution x = C/a 4 even for T = 0, and x can a °t 
as "mirage" radiation. The novel feature for T ^ is that this "mirage" effect appears 
through the decay of real brane matter, even if it was absent in the beginning (xi = 0)- 



3.3 The case of radiation for energy influx 

In the case of radiation the general solution of equations ( |3.l| )-( ^3|) was derived in the 
previous subsection and is given by equations ( glcj) , (Hull) . The expansion is that of a 
radiation-dominated universe with constant energy (xi I (27) + /?i) a i P er co-moving volume. 
The "mirage" energy density is diluted ~ a -4 . 

The explicit dependence on time will be discussed next in the case of flat space (k = 0), 
in which the energy density satisfies 



dp 2 
-dt + l P 



-T{p) • 



(3.15) 



If T(p) < for all p, and the "friction" term in the left hand side becomes suppressed for 
t — > 00, we expect an unbounded increase of p in this limit. For p > j//3 the low energy 
approximation employed in this section breaks down. The full treatment necessary in this 
case will be given in the next section. 

The actual situation is rather complicated and the details depend crucially on the form 
of T(p). Assuming that T(p) = Ap v with A < 0, the exact solution of equation ( 3.15 ) for 
v + 1,3/2 is 



I T?\ 1 — V 1 



2v 



I3-2i/ 



t 



). 



(3.16) 



where t = \A\t. For v = 1 the solution is 



P - - 

i „t-U 



Pi ¥ e 



and for v = 3/2 



-1/2 



(PiPi) 



-1/2 



2 V u 



(3.17) 



(3.18) 



Having neglected the term f3p from equations (|3.2[) and (3.3), these solutions are valid only if /3pf -C H, 



As we shall point out below this condition is eventually satisfied for all solutions in the case of outflow of 
the form discussed here and k = 0. 
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For < v < 1 we have p ~ for t — > do. For v = 1 the increase of the energy 

density at large t is exponential moderated by a power. For 1 < v < 3/2 the energy density 
diverges at a finite time 



r 2 *=r 2 *+^o^) 1 ~ 1 ' ■ (3-i9) 

A similar divergence appears for i/ = 3/2. For f > 3/2 a divergence occurs if the quantity 

D = {PM) 1 -" ~ ^-/r 3 (3-20) 

is negative. In the opposite case pP — > 1/D for i — ► oo, and the energy density diminishes: 
p ~ i~ 2 . 

As we discussed earlier, it is physically reasonable that the energy influx should stop at 
a certain value p cr , and be reversed for larger energy densities. The dynamical mechanism 
that localizes particles on the brane cannot operate for arbitrarily large energy densities. 
This modifies the solutions above that predict an unbounded increase of the energy density. 

A final observation that will be encountered again in the next section is that, despite 
of the fact that the energy density in most cases increases for large times, it can decrease 
at the initial stages. This is obvious from eq. ( |3.15 ). If at the time t r that the brane enters 



a radiation dominated era \T(p)\/p < 2/t, the energy density decreases for a certain time. 
3.4 Non-flat solutions 

In addition to the analytical special solutions discussed above, we would like to present a few 
suggestive numerical results concerning the k = ±1 cases. For u = 1, we substitute (|3.12j ), 
true for any A, into the second order equation ( |2.17| ), to obtain 

- + % + + + 2,w)(3p\ e - (1 - 3w)j Pl — ^ttt — c = . (3.21) 

a a 2 a 2 a 6(i+«0 > irL a 3(i+«0 ^ ' 

It is obvious from ( |3.2l| ) that for outflow with k = 1 and w > —1/3 the Universe will 
exhibit eternal deceleration. In particular, in the case of dust, figure |l| depicts the solution 
for a{t) of (^2i|) for some initial conditions for a, p. Notice that a(t) after a period of 
decrease, starts increasing again, thus delaying its re-collapse. Of course, with appropriate 
initial conditions one also obtains solutions with the standard FRW behavior. 

For an open (k = — 1) Universe with A > 0, w = 0, v = 1 we have found numerically a 
solution where p(t) monotonically decreases to zero, while a(t) starts with deceleration, but 
later on accelerates eternally; more specifically, q(t) — > + for t — * 0. Another possibility 
for k = — 1 and A < 0, w = 0, v = 1 allows for a Universe starting with acceleration at 
infinite densities, which later on, turns to deceleration with p approaching a constant value. 

4. General remarks on the outflow/influx equations 

In [ffl] the interested reader may find a rather detailed study of the fundamental equations 



of cosmology (2.19)-(2.21) in the particular context presented here. 
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Figure 2: Outflow, fc = 0, w = 0, v = 1. The arrows show the direction of increasing scale factor. 

The general analysis, however, was done numerically for various cases of the energy 
exchange functional dependence on the density p, various values of k and various types of 
matter on the brane. Two examples with of the resulting cosmologies are the ones shown 
on Figures 2 and 3. 

Figure 2 shows the global phase portrait of q = a/ a versus p for the case of dust on a 
flat (k = 0) brane, with energy outflow and T = Ap. The arrows point in the direction of 
increasing time, or equivalently of increasing scale factor for the expanding solutions we are 
interested in. In all cases the density decreases with time and the universe evolves towards 
the fixed point with vanishing energy density and zero q. One may, however, recognize 
two classes of evolutions. The one in which the universe decelerates at all times, and the 
one in which the universe first decelerates, then it lives a period of acceleration before it 
decelerates again to approach the fixed point form below the p axis. The behavior near the 
fixed point may be derived analytically. Also analytically one may derive the presence of 
the limiting parabola q = jp — Pp 2 , shown with a dotted line. For comparison, we also show 
the straight dashed line q = — (1 + 3w) r yp, which represents the standard FRW behavior. 

A lot richer and diverse is the case of energy influx. The cosmological evolution of 
dust on a flat brane with linear exchange function T{p) is presented in Figure 3. Again, 
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Figure 3: Influx, k = 0, w = 0, v = 1. 



as in the previous case, one may recognize the limiting parabola, restricting all possible 
histories of the brane. Its slope at the origin is also determined analytically and verified 
numerically. 

Not shown is the line q = — (1 + 3w) r yp of FRW cosmology. Due to the energy influx, 
the origin is not anymore a future fixed point. Instead, there are two fixed points, both 
with constant acceleration and constant energy density. The one on the left is an attractive 
fixed point for almost all initial conditions. The other is a saddle point corresponding to 
unstable trajectories. Clearly, there is now a variety of possible histories, but it is premature 
to concentrate on any one of them in an effort to understand our Universe. One has to 
include the dynamics of the bulk and solve the complete system, before such an attempt 
can be anything but academic. 

5. A first step towards a complete description of the brane- world cosmol- 



The analysis so far was based on the assumption that the bulk is some kind of reservoir 
and does not change appreciably during the cosmological evolution of the brane. So, the 
bulk was neglected and the autonomous system of equations for the brane and the energy 
exchange with the bulk was discussed. However, realistic applications of this scenario call 
for a more general treatment, including the dynamics of the bulk. 

A first step towards a complete analysis of the brane-bulk coupled system, is to try 
at least to describe the stable fixed point of the influx case, the far future behavior of the 
brane- world. 

It is straightforward to convince oneself that the metric 



satisfies the five-dimensional Einstein equations with energy-momentum tensor, whose non- 
vanishing components are (8ttG = 1) 



ogy. 



ds 2 = -dt 2 + e bt+c ^dx i dx i + dz 2 
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Tl = Tl = Tl = \{c 2 -l?) + 2c5(z) (5.3) 
T| = |( c 2- 26 2 ) (5.4) 

Ti = \bce{z) (5.5) 

To decide about the nature of matter in the bulk, one has to transform to a coordinate 
frame with vanishing Tg . This is achieved by a Lorentz transformation in the t—z subspace. 
Such a transformation will leave the —dt 2 + dz 2 part of the metric invariant, while the 
exponent in the metric will depend linearly on only one of the new coordinates. Specifically, 
one distinguishes two cases. 

Case A. |c| < The transformed metric takes the form 

with the corresponding energy-momentum tensor 

3 

T'(; = --(6 2 -c 2 )di a g(l,l,l,l,2). (5.7) 

One interpretation is that in the bulk there exists a perfect fluid with constant energy 
density and pressure, given by 

PB = ^(b 2 -c 2 )>0,P B = -\{b 2 - c 2 ) = -p B , P T = \(c 2 - b 2 ) = -2 PB (5.8) 

flowing towards the brane in the brane rest frame. 

Case B. |c| > \b\. The transformed metric takes the form 



ds 2 = -dt' 2 + e b ^ IZ ^ t 'dx i dx i + dz' 2 (5.6) 



ds 2 = -dt" 2 + e c ^ l ~ b2 l c2z "dx i dx i + dz" 2 (5.9) 
with the corresponding energy-momentum tensor 

T^ = -^6 2 -c 2 )diag(2,l,l,l,l). (5.10) 
interpreted as a perfect fluid with constant energy density and pressure, given by 

PB = 3 ~(b 2 -c 2 )<0,P B = -\(b 2 - c 2 ) = -\p B , Pt = ~\{b 2 - c 2 ) = -\p B (5.11) 

again flowing towards the brane in the brane rest frame. 
The induced metric on the brane is just 

ds 2 = -dt 2 + e^dx^x 1 , (5.12) 

which for b > is an exponentially expanding deSitter space-time with constant q. 

In addition, on the brane there is a perfect fluid again with constant energy density 
and pressure 

p b = -3c ,p b = 2c = -2p b /3. (5.13) 



-13- 
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For c < 0, to avoid exponential blow-up of the metric at transverse-space infinity, the 
brane density is positive. Assuming linear state equation for the matter on the brane, the 
situation corresponds to Wf, = —2/3. 

Finally, there is constant energy influx from the bulk. Indeed, T{z > 0) = T®(z > 
0)/2 = 3 be/ A < 0, for the chosen signs of b and c. 

It seems that the above metric provides a complete description of the brane-world 
cosmology at late times. It realises the main characteristics of the stable fixed point of 
Figure 3. It describes a brane carrying constant positive density, expanding with constant 
acceleration in the presence of constant influx from the bulk. 

There are however a few unsatisfactory features in the above solution. First, instead 
of dust on the brane we had to introduce a fluid with w = —2/3. Second, assuming a linear 
equation of state for the matter in the bulk, the above solution of 5-dimensional Einstein 
gravity with a brane does not satisfy in the bulk any of the known energy conditions. 
Indeed, the perfect fluid in the bulk does not satisfy the null energy condition (NEC) 
(pb + Pb > and pb + Pt > 0) and consequently, it does not satisfy any of the other known 
energy conditions (weak, strong, dominant or null dominant). Nevertheless, it should be 
pointed out in connection with these comments, that on the one hand violation of the 
energy conditions does not necessarily mean that the theory is fundamentally sick and 
unacceptable, by being acausal or unstable (Remember, for instance, the recent discussion 
in the literature of the so called "phantom" matter). On the other hand one does not have 
to interpret the above solution assuming a linear relation between pressure and density (see 
for instance |21[|) and this might allow for an interpretation with more "natural" matter, 
both on the brane and in the bulk. 

6. Conclusions 

A rather detailed mathematical analysis of the role of the brane-bulk energy exchange 
on the evolution of a Brane Universe, adequate for a wide range of potentially realistic 
implementations, was presented. The effective brane cosmological equations were derived 
with perfect fluid matter on the brane, constant energy-momentum tensor in the bulk 
and non-vanishing exchange between them. A detailed study of the solutions of these 
unconventional equations was performed in the case of zero effective cosmological constant 
on the brane, in order to reduce to the Randall-Sundrum vacuum in the absence of matter. 
The analysis revealed a rather rich variety of possible cosmologies, depending upon the 
precise form of the exchange term, the topology of 3-space, and the nature of matter on 
the brane. 

General characteristic properties of the solutions were presented, together with a few 
special but particularly interesting solutions, obtained in the limit of low energy density 
on the brane. One of these is the exactly solvable case of radiation, where a mirage ra- 
diation effect appears through the decay of real brane matter. Another, is a De Sitter 
fixed point solution, obtained in the case of energy influx, even without pure vacuum 
energy, which in addition is stable for a wide range of reasonable forms of energy ex- 
change. 
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Clearly, the formulation of a detailed brane-world cosmological scenario requires a 
complete analysis of the full brane-bulk system. A first step towards understanding the 
full dynamics was taken above and a description of the asymptotic times cosmology was 
given. It has a positive brane density and a positive q, and represents nicely the attactive 
fixed point plotted in Figure 3, for the case of energy influx. Many open questions, such as 
an all-time description of the brane-world evolution, the duration of accelerating periods, 
the creation of primordial fluctuations, and the compatibility with conventional cosmology 
at low energy densities, should be addressed. However, we believe that the cosmological 
evolution in the context presented here has many novel features, that may provide answers 
to outstanding questions of modern cosmology. 
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